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Abstract 

In order to explain the non-hierarchical neutrino mixing angles and the absence of 
lepton flavor violating processes in the context of warped extra dimensions one needs 
to introduce bulk flavor symmetries. We present a simple model of lepton masses and 
mixings in RS models based on the A4 non-abelian discrete symmetry. The virtues 
of this choice are: (i) the natural appearance of the tri-bimaximal mixing pattern; 
{ii) the complete absence of tree-level flavor violations in the neutral sector; {Hi) 
the absence of flavor gauge bosons; (iv) the hierarchies in the charged lepton masses 
are explained via wave-function overlaps. We present the minimal field content and 
symmetry breaking pattern necessary to obtain a successful model of this type. The 
bounds from electroweak precision measurements allow the KK mass scale to be as low 
as ~ 3 TeV. Tree-level lepton flavor violation is absent in this model, while the loop 
induced /x — > branching fraction is safely below the experimental bound. 



1 Introduction 



Warped extra dimensions [1] provide a simple framework for fermion masses: exponential 
hierarchies are naturally generated due to the overlap of fermion and Higgs wave functions [2, 
3], implementing the split fermion idea of [4]. In the simplest 5D "anarchic" approach [5,6], 
where both the 5D bulk masses and the brane Yukawa couplings are assumed to be random, 
one generates a hierarchy both in the 4D standard model (SM) fermion masses and their 
mixing angles. This seems to fit the observed pattern of quark mixings and masses very 
well, since both the CKM matrix and quark masses show a hierarchical pattern. The lepton 
sector, however, is different: two of the observed neutrino mixing angles are close to maximal, 
rather than being hierarchical [7]. This suggests that one needs to radically change the 
implementation of fermion masses in warped extra dimensional models in order to achieve 
the correct lepton mixing pattern. Instead of a fully 5D anarchical approach, it calls for 
partial flavor symmetries, which will make sure that two of the neutrino mixing angles are 
not small, but close to maximal. The usual hierarchies can still be used to generate the 
hierarchies in charged lepton sector. 

In fact, the necessary appearance of flavor symmetries is welcome for these models. An- 
archic 5D flavor structure necessarily gives rise to flavor changing neutral currents (FCNC's) 
in the quark sector, and to lepton flavor violation (LFV) in the lepton sector at tree- 
level [3,5,6,8]. While these flavor violations are suppressed by the so-called RS-GIM mech- 
anism [6,8], the KK mass scale still has to be as large as 20 TeV in order to safely suppress 
FCNC's in the quark sector [9]. In the lepton sector the anarchic approach [10] imposes a 
bound of order 10 TeV on the KK mass scale [11]. These bounds imply that the theory is 
out of reach for the LHC and unlikely to be useful for solving the (little) hierarchy problem. 
Reduction of these bounds require additional flavor symmetries [12,13]. For the leptons, we 
have seen that such flavor symmetries are necessary to begin with, to protect the mixing 
angles from becoming hierarchical. The aim of this paper is to show that we can indeed use 
these flavor symmetries in the lepton sector to eliminate the LFV bounds on the model, and 
at the same time get a realistic pattern of masses and mixings. Recently, an alternative ap- 
proach to lepton masses in RS was proposed in [14] (following the suggestion of [15]), which 
however does not readily explain the appearance of the large (non-hierarchical) neutrino 
mixing angles. 

For our model, we pick the most popular global symmetry used in neutrino model build- 
ing, an A 4 discrete symmetry [16, 17]. This, by the virtue of being discrete, has the added 
benefit that no additional gauge bosons have to be introduced, even if the symmetry is 
gauged. From the many studies of the implications of the A 4 symmetry in 4D models we 
expect that this symmetry can give the correct mixing structure. The mass hierarchies in 
the charged sector can still be generated as usual in RS models via fermion overlaps. The 
mass hierarchy in the neutrino sector is not that large, and can be readily incorporated by 
choosing 0(1) factors in the neutrino mass matrix. Finally, since A 4 is a non-abelian discrete 
symmetry it has the potential of greatly reducing bounds from LFV. The reason for that is 
that by using a non-singlet representation under A 4 we can force the bulk wave functions of 
some of the left handed fields to be universal. 



1 



Indeed, we find that with appropriate choice of representations, the tri-bimaximal mixing 
pattern [18] characteristic of the A4 symmetry can be reproduced by the leading order 
operators. Higher order terms can result in non-zero #13, while maintaining the predictions 
for 8 12 , 6*23 within the experimentally allowed range. Most importantly, our A 4 based model 
eliminates all tree-level sources for LFV. In fact, LFV is completely absent in the neutral 
current interactions, and shows up only through charged currents involving neutrino mixing. 
In the original RS neutrino mass model of [2] these loops turn out to be problematic [19], 
since the interaction of the heavy KK neutrinos with the SM fields is unsuppressed. In our 
model, by putting the charged leptons in the bulk (and peaked on the UV brane) these 
couplings are strongly suppressed and the bounds from loop induced decays are significantly 
reduced. The most important experimental bounds in the lepton sector are those from the 
electroweak precision (EWP) constraints. They are, however, quite mild, with KK masses 
of order 3 TeV generically allowed as in most other RS models. 

The above arguments show that introducing the discrete non-abelian lepton flavor sym- 
metry greatly improves over the generic RS lepton flavor models. Moreover, they are also 
improving the straight 4D implementations of A 4 neutrino models in several aspects. First, 
they explain the hierarchy in the charged lepton sector. Second, by putting one A4 breaking 
VEV on the UV brane (breaking the group to Z2), and the other on the IR (breaking it to 
Z3) questions regarding vacuum alignment are eliminated. Finally, the appearance of the 
correct right handed neutrino mass scale (somewhat below Mpi and Mqut) can be explained 
by the partial compositeness of the right handed neutrino. 

The paper is organized as follows: in Section 2 we give the general setup, introduce the 
A4 representations and calculate the mixing matrices at leading order. In Section 3 we show 
the effects of higher dimensional operators on the mixing angles. In Section 4 we present 
a numerical scan of the parameter space and discuss the electroweak precision bounds. In 
Section 5 we show that LFV is completely absent at the tree-level in this model, and estimate 
the loop induced \x — > ej branching fraction. We conclude in Section 6. 



with a UV brane at z = R (R is also the AdS curvature scale) and an IR brane at z = Rl . 
The magnitude of the scales is given by R~ l ~ Mpi and R'~ l ~ 1 TeV. The electroweak 
gauge group is extended to an SU(2)lxSU(2)#xU(1)b-l gauge symmetry in the bulk to 
incorporate custodial symmetry [20] . This symmetry is reduced on the UV brane to the SM 
group SU(2) £ xU(l)y, while it breaks down to SU(2) Z) xU(1) j b-l on the IR brane. 

The matter content is summarized in Table [TJ We assume that there is a separate doublet 
for every SM lepton, including the right handed neutrino: an SU(2)^ doublet, for every 
left handed doublet L, a separate SU(2)« doublet ^ e ,n,T for every right handed charged 



2 The Setup 



We are assuming an AdSs bulk metric 
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Table 1: Fields and their gauge and flavor charges. 



lepton, e, /i, r, and a right handed doublet ty v for every right handed neutrino v. A 5D 
fermion correspond to 2 Weyl fermions of opposite chirality in 4D: 

•-(*)• (2 - 2) 

In the KK decompostion, 4D chirality follows from the boundary conditions at the two 
end points of the extra dimension. The Lorentz structure forces x an d ip to have opposite 
boundary conditions while a massless mode appears only for Neumann boundary condition at 
both ends. For a complete description of fermionic boundary conditions see [21]. The profile 
of the would-be zero mode is then entirely dictated by the 5D mass of Conventionally, 
this mass is normalized as c/R and for c > 1/2 (resp. c < —1/2), a x _zero mode (resp. 
ip-zero mode) is exponentially localized on the UV brane. In our setup, we assume that the 
boundary conditions (in the absence of the localized mass terms) are chosen as: 

v L = ( l l+ , +l ) » w = ( j+. :])». = (£[;;:]) (2.3) 

where [±] refers to a Neumann (Dirichlet) boundary condition on both branes for the \ 
component, while the ip ones simply have the opposite conditions. Hence there is a left 
handed zero mode for the left handed doublets in and a single right handed zero mode 
in ty e ,ti,r and ty u eachQ These fields have bulk masses (in units of the AdS curvature) given by 
Cl, c e and c u . For the general case these would be hermitian 3x3 flavor matrices. This is not 
the case here as we impose an A4XZ2 global symmetry in the bulk of the theory. In order to 
get the correct neutrino mass spectrum, we assign the three charged lepton doublets as well 
as the three right handed neutrinos to the 3 dimensional representation of A4. Thus for these 
fields there is just one common c-parameter each: ci and c v . On the other hand, we assign 
the right handed charged leptons to the three inequivalent one dimensional representations 
of A 4 : l+l'+l". Thus there are three separate c's in this sector: c e , c M and c T . The purpose 
of the Z2 symmetry is to eliminate certain brane localized operators that would otherwise 
contribute to the mass matrices at leading order. 

1 It is well known that the [H — ] boundary conditions (for a x-type Weyl fermion) can lead to an extremely 
light KK-state for c > 1/2. Here we are safe from this problem as the right handed zero modes which satisfy 
this type of boundary conditions are localized close to the UV brane, i.e. c < —1/2. 



3 



The symmetry breaking is achieved via brane localized scalars. Besides the SM Higgs, 

H , that is localized on the IR brane, we assume the following scalars to break the discrete 
symmetries: <fi on the UV brane and <$' on the IR brane, both of which are in the 3 of A4. 
We assume that these two scalars develop VEVs in different directions: (ft breaks A 4 to Z 2 , 
while cf)' to Z 3 . This is achieved by the VEVs 

(<f>') = (v , ,v , ,v'), = 0,0), (2.4) 

in the basis where the generator corresponding to generator S of A4 is diagonal (see Ap- 
pendix A for summary on A 4 ). Note, that once such a basis is chosen, these are the most 
general VEVs which preserve Z3 and Z2 subgroups of A4 respectively, up to a trivial permu- 
tation of the basis. 

We now write the most general Yukawa terms respecting both gauge and flavor symme- 
tries (in addition one needs to write localized kinetic and potential terms for the localized 
scalars) : 

Cuv = _ ^VvVv - + h.c. H , (2.5) 

C IR = -p L H* u - * (*,0') H* e - ^ PM')' H% - % pM')" H* T + h.c. + ■ ■ ■ , 

where \I/ T = (x, V>) and ■ ■ ■ stands for higher dimensional operators. We use the notation 
of [17] for writing A^invariants: (),()' and ()" denote the terms of 3 x 3 that transform as 

I, 1' and 1" respectively, see appendix |A} 

Once the electroweak and A 4 symmetries are spontaneously broken, the boundary terms 
lead to boundary conditions mixing all fermions. Then, the spectrum is obtained by solving 
the bulk equation of motion in the presence of these mixed boundary conditions. The light 
modes, however, are quite insensitive to the boundary terms, and so they can be treated 
as a small perturbation. Hence to leading order the low energy spectrum may be obtained 
by using the zero mode wave functions. This defines the Zero Mode Approximation (ZMA) 
whose accuracy depends on how light the light masses are. As the largest mass is that of 
the r, about 1 GeV, the ZMA turns out to be as accurate as m T R' ~ 10~ 3 for zero mode 
masses. 

To follow the conventional RS literature we introduce the RS flavor functions / and F; 
these give the wave functions of the zero mode fermions on the IR and UV branes: 



fR>\RJ \R'J VR\R 
with 

f = VrrTc F = 

y/l ~ (R/R') 1 - 2 ^ C y/l ~ (R/R'yc- 1 ' 



The IR boundary terms of (2.6) lead to the following Dirac mass matrices for the zero mode 
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charged leptons and neutrinos: 




V2RA \Vef-e " 2 yJ-» 

M » = vJlU vMI' 1 ° 1 0| ' (2 ' 9) 

where we have introduced the shorthanded notation fa = f c . and f-i_— f- Ci ), w is the cubic 




root of unity, lo = e 2m ^ 3 , and vh ~ 250 GeV. The UV terms of (2.6) generate a Majorana 
mass matrix for right handed neutrinos of the form: 

(e s 0\ 

M M = F 2 _ U R- 1 e s eA. (2.10) 
\0 e t e s J 

with e s = M/A and e t = x u v/A. 

Note, that while the mass hierarchy in the charged lepton sector is generated via the wave 
function overlaps in the usual way, the non-degeneracy in the neutrino sector is actually due 
to the two different kind of Majorana term allowed on the UV brane. As discuss below, 
the required neutrino mass hierarchy will actually require e s ~ e t . This is more naturally 
achieved if the singlet Majorana mass is actually also originating from an operator with a 
singlet scalar field VEV. This could for example be enforced by imposing an additional Z3 
global symmetry, and an additional scalar field £ with VEV (£) = M. 

After integrating out the heavy right handed neutrinos, one ends up with a seesaw type 
Majorana mass matrix of the left handed neutrinos 

Mm = -M» D ■ (M^y 1 ■ (MQ T 
'2A' 2 R' 2 F\~ 



-vl ZfJ L l- U [ e s /A -e t /A\, (2.11) 




with A = e 2 — e 2 . The diagonalization procedure is the same as that of usual A 4 four- 
dimensional models. The charged lepton mass matrix becomes diagonal once the left-handed 
fields have been rotated as L — ► V • L with 



V = —= 1 u 2 u , (2.12) 




a symmetric unitary matrix. We emphasize already here the most important properties of 
the A 4 mass matrices: the left handed rotation needed for diagonalizing the mass matrix 
is independent of the actual magnitudes of the masses, and the right handed rotation is 
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the unit matrix (i.e. no right handed rotation is necessary). Therefore, the charged lepton 
masses are 

'Vef-e 



\*-M 



D 



h 



V2R'A' 2 








Unf—n 






Vrf-7 



(2.13) 



The charged lepton mass hierarchies are due to the hierarchies on f- e ,-n ,-r and the A 4 
embedding of the right handed charged leptons allows for three different c's which can be 
set to generate the physical charged lepton masses. 

Next we move to the light neutrino sector. We work in the basis of diagonal charged 
leptons which is obtained by performing the rotation on the entire left-handed doublet with 
V. Then, the light neutrino Majorana mass matrix is diagonalized by the Harrison-Perkins- 
Scott (HPS) matrix [18] 

/ v^Ts 1A/3 \ 

U H ps= -l/Ve 1A/3 -1/V2 , (2.14) 
\-l/VQ 1/V3 1/V2 J 

corresponding to 9±3 = 0, sin 2 (2# 12 ) =8/9 and 9 23 = ir/4. This tri-bimaximal mixing matrix 
is very close to the the best fit obtained from present oscillation data. The neutrino mass 
eigenstates are 

i 



U 



HPS 



• V* • ■ V* • U HPS = 
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(2.15) 



where the overall mass scale is set by the combination 



m 



'2A' 2 R' 2 F 

The neutrino mass-squared splittings are given by 
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with r = e t /e s 
see that r solves the following algebraic equation: 

x — 1 



(1-rf. 

Combining the last two relations we 
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Am 2 sol /Am 2 atm for |r| < 2 or x 
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(2.19) 



—Am 2 sol /Am 2 atm when |r| > 2. Finally, e s is 



found by inverting one of the two relations, for instance, 
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Table 2: Approximate numerical values of the neutrino masses and UV VEVs for the 4 
possible solutions of r = e t /e s . The masses are given in units of 10~ 3 eV. 



When we impose the measured values [7] for the mass splittings 

Am 2 sol ~ 7.9 x 10" 5 eV 2 , Am 2 tm ~ 2.6 x 10~ 3 eV 2 , (2.21) 

we find four solutions for the neutrino mass spectrum corresponding to 

r « {0.79, 1.19, -2.01, -1.99}. (2.22) 

The first (last) 2 solutions lead to a normal (inverted) hierarchy spectrum. Another general 
feature of A4 is a prediction of the absolute neutrino mass scale, m. We find that the mass 
of the heaviest neutrino ends up being slightly above the atmospheric splitting. 

We close this section by presenting a set of numerical values for the Lagrangian parameters 
which reproduce the mass spectra. The brane positions are R^ 1 = 10 19 GeV and R'^ 1 = 
1.5 TeV, in order to keep the KK gauge bosons (with itlkk = 3 — 4 TeV) in the reach of the 
LHC. The Higgs VEV turns out to be vh = 255.5 GeV, which is obtained after matching 
the bulk gauge couplings such that the weak boson masses as well as the fine-structure 
constant at the Z pole take their physical values: m,\y = 80.40 GeV, mz = 91.19 GeV and 
a em( m z) — 128. As usual [22], we get a Higgs VEV slightly above the SM value due to the 
suppression of the W, Z wave functions on the IR brane and the additional contributions to 
the gauge boson masses from the wave function curvature terms. 

The charged lepton masses are reproduced using the following choice of parameters: 
ci = 0.51, c e = —0.75, c M = —0.59, c T = —0.51, y e = 1.53, y^ = 1.55, y T = 3.04. Indeed, 
together with A' = R 1 ' 1 and v'R' = 0.1, one gets m e = 0.511 MeV, m M = 106 MeV and 
m T = 1.77 GeV. We show later on that the model defined with this set of parameters passes 
all leptonic electroweak precision bounds while perturbation theory remains under control 
up to E = 3m,KK- As we saw above, in the neutrino sector the solar and atmospheric 



mass splittings, Eq. (2.21), reduce the ratio r to only four possible values, Eq. (2.22). The 
corresponding mass spectra, as well as the Majorana masses on the UV brane (for A = 
x u = 1 and c v = —0.37), are reported in Table |2] 



3 Higher Order Corrections and Perturbativity Bounds 

After spontaneous breaking of A 4 the lowest dimensional boundary terms of ( |2.6[ ) generate 
the tri-bimaximal pattern for the mixing angles as well as the charged lepton and neutrino 
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mass hierarchies. In order to consider this construction as really meaningful, it is necessary 
to study its stability under corrections from higher dimensional terms on the branes as well 
as radiative corrections. Another motivation for looking at deviations from tri-bimaximal 
mixings is the possibility to get a non-zero 6^3 , in case it turns out to be non- vanishing 
experimentally. 

We focus first on the UV brane and show that the effects of higher dimensional A 4 
invariants lead to the same pattern for the Majorana mass matrix, except that some entries 
become complex. This is the only source of non-zero #13 in this model. We start with writing 
down the most general higher order operators allowed on the UV boundary: 

An 

- 5C UV = ^2\ n —ij u ij u + h.c. (3.1) 



n>2 



with n insertions of 0. Now due to the Z 2 -preserving VEV ((f)) = (v, 0, 0), it is straightforward 
to show that 3 transforms as under A 4 . Thus all the higher order effects that cannot be 
reabsorbed into a redefinition of the lowest order parameters arise from one operator: 

<f) 2 

- SCuv = A 2 -^M^ + h.c. (3.2) 
This term contains actually three independent A 4 invariants which lead to complex diagonal 



elements of the Majorana mass matrix and the pattern of (2.10) is to be replaced by: 



(e s + 5 1 
M U M = F 2 u R~ l e s + 6 2 e t ] . (3.3) 
V e t e s + 8l 

with 5i ~ 0(v 2 /A 2 ) and 5\ real. The complex entries induce both a deviation of #i 2 from its 
maximal value and a non-zero sin(#i3) of 0(5). 

On the IR brane, we show that the higher order corrections can only affect 6*12, as long 
as Z 3 remains unbroken on this boundary. Again we start with the most general higher 
dimensional invariants which take the following form: 

- 8C IR = E E Kn*^H% + £ ^l^H^ + h.c. (3.4) 

Since (<f)') = (v',v',v') is Z3 symmetric, <f)' 2 transforms as l + <p'. Then non-trivial corrections 
reduce to terms with only one <// insertion: 

- 5C IR = ^ L 4tH* v + h.c. (3.5) 

Since it is suppressed by only one power of (v'A') compared to the lowest order terms, this 
operator may easily destabilize the HPS pattern. The extra Z 2 flavor symmetry is useful 
here, as we now discuss. We choose the A 4 triplets odd under this additional Z 2 , and then 
this operator is forbidden by the Z 2 symmetry. However, since <fi' 2 = 1 + 0', the next Z 2 
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even operator generates the same type of correction as the one linear in 0' but with a higher 



suppression factor. In that case (3.5) has to be replaced by 



- 5C IR = ^ L <p' 2 HV u . (3.6) 

This operator corrects the Dirac neutrino mass matrix by introducing off diagonal elements 
of 0[v' 2 /A /2 ). As this term contains three independent A4XZ2 invariants, the diagonal Dirac 
matrix of ( 2.9[ ) becomes: 




•Md = Vvhf- V j=^ I e 3 l + e 1 6 2 |, (3.7) 

with e.j ~ 0(v' 2 / A' 2 ). One can easily check (see Appendix B) that only sin(6 , 12 ) is affected 
and its deviation from the HPS prediction is of 0{v' 2 / A' 2 ). 

There is another important feature we want to stress at that point, which is the fact 
that there is no correction to the charged lepton mass matrix. This means that the ^ e ,fi,r 
fields need not be rotated again to get to the diagonal charged leptons basis even in the 
presence of these higher order terms. As we shall see in section 5, an immediate consequence 
is the absence of tree level LVF in this model, even when higher dimensional operators on 
are considered. 

We close this section by presenting an NDA estimates for the allowed sizes of the IR brane 
localized operators. These bounds are important for estimating how large the deviations from 
the HPS mixing matrix could actually be. For this purpose, we write down again the most 
general Lagrangian on the IR brane 

- C IR = V ^ L H^ U + + ^^ L H<j)' 2 ^ u + h.c. (3.8) 

We require that the theory remains perturbative up to scale En = Nttikk, which corresponds 
to the first A^ th KK modes of the theory being weakly coupled. We should require N > 3, 
such that at least the first few KK modes can be treated in perturbation theory. We can 
then systematically require that by the time this energy En is reached all loop corrections 
are still smaller than the lowest order terms. For example, for the first operator above there 
is a one loop correction to the Yukawa vertex itself (Fig. [l^,), which implies 

Vl E 3 N ^ , E 



< yu^, (3.9) 



16tt 2 A' 3 - *" A 

where the linear running of the coupling has been taken into account. Recalling wikkR' — 2, 
this implies for A' = R'~ l the perturbativity constraint 

Vu<2. (3.10) 

Similarly, for the second operator there is a two-loop correction to the tree-level operator 
(Fig. lib), we get 

Ve N <vM - y e <4 (3.11) 



(1 6vr 2)2 A /6 - ^ A /2 
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(a) (b) (c) 

Figure 1: Loop diagrams correcting the operators localized on the IR brane which are 
relevant for NDA estimates. 

again for A' = R'~ l . Finally, the third operator gives a one-loop correction to the first 
operator (Fig. fit), which implies the relation 



K 2 E% E N . . 

16^ A*- Vv A 7 " ~* ^<4y,<8. (3.12) 



We can then estimate that the higher dimensional terms corrections to sin(# 12 ) are suppressed 
compared to the lowest order term by at least A{y'R') 2 . 

We now apply the same arguments on the UV brane where we have 

- £>uv = |t"^Vv + Xvt^v^v + A 2 ^?/Wv + h.c. (3.13) 

We require that the theory remains perturbative up to the natural scale on that boundary, 
namely until E = R" 1 . First of all we focus on the last operator which contribute at one 
loop to its own vertex (Fig. |2]i) which implies the usual constraint 

A 2 <4vr, (3.14) 

where we assumed A = R^ 1 . On the other hand, the second operator induces one-loop 
corrections to the mass M, the third operator, as well as its own vertex. From the diagrams 
of Figs. [2£i, [2]d and ^ we derive the following relations: 



x v < Atc^ s , x v < Xl fi V^ < (4vr) 3/4 ~ 7, x v < 4tt. (3.15) 



We showed in the previous section that in order to reproduce the observed hierarchical 
neutrino mass splittings, one must have e s ~ e t with e t = x u (vR). Hence assuming a not 
so small suppression factor vR ~ 0.1, such that #13 is not dramatically tiny, the first of the 
above relations rewrites as x v < 16ir 2 (vR) ~ 15. Thus we end up with the perturbative 
constraint 

x v < 7. (3.16) 
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/ \ / \ 

/ x ^ ^ x \ 



(d) 

Figure 2: Loop diagrams illustrating the NDA estimates on the UV brane. 

4 Numerical Scans and Electroweak Precision Bounds 

One of the main constraints in models with new physics at the TeV scale are the electroweak 
precision measurements (EWPM). For RS models with custodial symmetry the generic bound 
on the KK scale is about tukk > 3 TeV, mostly from the contribution to the S-parameter [20] . 
Here we will check that the lepton sector of our model indeed passes these tests for KK scales 
of order 3 TeV. El 

The simplest way of checking the electroweak precision constraints in an RS model with 
bulk fields is to canonically normalize the SM gauge fields, and to determine the parameters 
g 5 , g' 5 and vh by requiring that the measured values of M W ,M Z and the electromagnetic 
coupling e are reproduced. This choice is somewhat unconventional, since M w is less accu- 
rately measured than Gf, however it simplifies the matching of the 5D parameters to the 
observables significantly. In this scheme all corrections to electroweak precision observables 
will be contained in the fermion-gauge boson vertices, which can be simply calculated by 
wave function overlaps, and compared to the SM predictions in terms of the above input 

1 It is possible to cancel the S-parameter by tuning all LH bulk fermion masses to be ~ 0.5 [23], as it is 
necessary in higgsless models [24]. In this case however one does not get any mass hierarchies and those 
need to be introduced by hand as in [12]. 
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Figure 3: Deviation from SM Z couplings of lz and Ir as function of c's. di = (g z — 
9z,sm) / 9z,sm and d R = (9z-9z,sm)/9z,sm are plotted in per mil units. We took R'~ l = 1.5 
TeV and used m w = 80.403 GeV, m z = 91.1876 GeV and e = e(/i = m z ) = ^4^jY2& as 
physical input observables. The red regions are excluded by EWPM in the leptonic sector. 



parameters, which for the lepton- Z-couplings are 



9z,sm 



9z,SM 



rn 



w 



rrir 



m z 



m z 

j 

mw 



m 



m w Jl - 



m w 
m| 



(4.1) 
(4.2) 



In a warped extra dimension the couplings of the left and right handed charged leptons are: 



9z - 



R,i 

9z 



1 



h" 



dz 



R' 



dz 



R 



R 



R 



[95La L ' 3 (z,m z ) + g 5 a x (z,m z )] Xc L (z)' 
[95Ra R ' 3 (z,m z ) +g b a x (z } m z )] ip^zf 



(4.3) 
(4.4) 



with a(z, m z ) the wave functions of the three neutral gauge fields containing the physical Z, 
while the various g§ s denote the bulk gauge coupling^] The origin of the deviations from 
the SM predictions is the non-flatness of the Z wave function. If the Z boson was massless, 
its wave function would be flat and the lepton couplings would become universal thanks to 
the orthonormality of their wave functions. However as soon as the A4 symmetry is imposed, 
the left handed lepton couplings remain flavor blind. In the ZMA the Z coupling depends 
only on the bulk mass parameter c^. As usual we assume that the SM leptons are localized 
on the UV brane, cl > l/2,c e;M>T < —1/2. Since the r is the heaviest lepton, it has to be 



2 For completeness we review in appendix [c] how this quantities are calculated in terms of physical ob- 
servables. 
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localized closest to the IR brane, so it will be most sensitive to the non-flatness of the Z 
close to the IR brane, and so one expect the biggest deviations in the r couplings. 

The couplings of the charged leptons have been measured very precisely at the LEP 
experiment [7]. Here we will require that all lepton couplings are within 0.2% of the SM 
prediction. The plots of figure [3] show the deviation of the charged lepton couplings to the 
Z from their SM values as a function of the c's. We see that, while the Cg »| Cciii be? clS close 
to —1/2 as required to reproduce the charged fermion hierarchy, cl cannot depart too much 
from 1/2 to remain within the experimental bound. The fact that cl is preferred to be close 
to 1/2 may be surprising at first, but we remind the reader that these vertex corrections 
also include the S-parameter. This is actually a welcome fact, since in order to keep the r 
Yukawa coupling perturbative we have to take cl close to 1/2 anyway. Thus we conclude 
from Fig. || that with a 3 TeV KK mass scale the electroweak precision constraints in the 
lepton sector are safely satisfied. 

Next we want to scan over the model's parameters and show explicitly that the neutrino 
mixing angles, which deviate from the HPS pattern in the presence of higher order operators, 
are within the allowed range of the existing results of the neutrino experiments. Once the 
usual RS solution to the hierarchy problem is imposed, we still have 12 free parameters in 
our setup: 5 c's, 5 Yukawas and 2 VEVs relevant for lepton masses v, v', while adding higher 
order operators brings 6 more coupling constants. We use the measured lepton masses and 
the best fit neutrino mass splittings to fix 5 of the lowest order parameters, which leaves still 
a lot of freedom to explore. Therefore we add some mild assumptions in order to simplify 
the parameter space and to try to only focus on the main predictions without having to go 
into the details of the structure of the higher dimensional operators. First of all we impose 
c T = —cl and keep cl as a free parameter. We also take the Yukawas on the IR brane 
to saturate the perturbative bounds: y e ^ tT = 2y u = 4. Thus imposing m T = 1.77 GeV, 
m e = 0.511 MeV and m M = 106 MeV in turn fixes v' and c 6)A , as functions of Cl- Hence all 
IR brane effects will be encoded in one parameter cl, which is constrained by the EWPM as 
discussed above. On the UV brane the solar and atmospheric splittings fix the ratio M/v and 
the overall neutrino mass scale. Then taking x v = 1 we end up with v = v(c u ) as the only free 
parameter on this boundary. Furthermore, (vR) 2 controls the size of #13 generated through 
higher order corrections. Finally we present in figure [4] some contours of the mixing angles in 
the plane (cl, c u ) when one among the possible combinations of higher order invariants are 
included on both branes. We considered separately the effect on 9i2 from the IR brane and 
the predictions for all angles from the UV brane. In order to demonstrate the robustness 
of the tri-bimaximal pattern under higher order corrections in this model we have selected 
the worst case where these operators saturate their perturbative limits, namely for A2 = 47r 
and k 2 = 8. We conclude from Fig. [4] that, even when the deviations from HPS angles are 
the largest possible, there is still a viable region satisfying the constraints sin 2 (2^13) < 0.19 
(90CL), sin 2 (2# 23 ) > 0.92 (90CL) and 0.73 < sin 2 (2# 12 ) < 0.95 (3a) [7,25]. Obviously the 
smaller the higher order terms are the closer one would get to the HPS pattern. 
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Figure 4: Scan of the parameter space reduced to (cl,c u ) as motivated in the text. The 
red regions are excluded by electroweak precision constraint on the Z coupling. We then 
show within this region some contours of the mixing angles delineating the largest values 
we typically obtain in the presence of higher dimensional corrections. The two horizontal 
contours are for the following values of sin 2 (20 12 ) = 0.90,0.95 (from bottom to top) where 
only the IR higher dimensional operator is added with e\ = e 2 = and e 3 = 8(v'R') 2 . 
The oblique lines are contours of sin 2 (26 , 13 ) = 0.01,0.19 (from left to right) and contours of 
sin 2 (2#i2) = 0.90, 0.95 (from left to right) generated by the higher order Majorana mass on 
the UV brane, for 81 = 5 2 = and <5 3 = 4tc(vR) 2 . 

5 Constraints from Lepton Flavor Violation 

Flavor models usually predict new sources of flavor violations, and so are only viable when 
the flavor scale is pushed to very high values. The flavor scale for quarks in the usual 
anarchic RS flavor models is around 20 TeV [9]. Thus it is very important to understand 
what the possible sources of lepton flavor violation (LFV) are in this model of lepton masses. 
Generically there are two types of LFV sources: tree-level LFV via the exchange of heavy 
neutral particles (like Z') or off-diagonal Z couplings, or loop induced rare decays via charged 
current interactions. 

Tree-level LFV operators have been considered in [11], and it was found that the lepton 
flavor scale is at least 5 TeV, and larger for smaller brane Yukawa couplings. We will first 
show here that the structure of the A 4 symmetry used here is such that all tree-level sources 
of lepton flavor violation are absent in this theory. This can be seen in the following way. By 
the choice of the A4 representation the wave functions of the left handed SM fermions are 
flavor universal (since they have the same Ci). So the only source of flavor violations in the 
charged lepton sector is the choice of different c ej/JjT necessary for the mass hierarchies, and 
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the brane Yukawa matrix (2.8). The couplings to the neutral bulk fields like the KK tower 
of the Z are controlled by the c ejAtjT , and will be non- universal. Thus any RH rotation to 
diagonalize the brane Yukawa matrix would induce tree-level LFV's. However, we have seen 
that one of the magic properties of the A 4 models is that the charged lepton mass matrix is 
diagonalized by a single left handed rotation, and no right handed rotation is necessary to 
diagonalize the mass matrix, as described in (2.13). The left-handed rotation is harmless, 
since the bulk wave functions are universal in the LH sector, while in the dangerous RH 
sector there is no rotation necessary. This implies that there is a basis where the kinetic 
terms and the mass terms for the charged leptons are simultaneously diagonal, and so as a 
consequence there is no tree-level LFV in this model. The lepton embedding in A 4 provided 
us with the necessary conditions to ensure the absence of LFV, namely universal left handed 
c's and the absence of redefinition of the right handed fields. Moreover we want to stress 
that this result remains unchanged once higher dimensional brane operators are considered 
as they were shown not to affect the rotation matrices of the charged leptons. 

Thus all lepton flavor violation arises from charged current interactions. In the SM loops 
involving neutrinos give extremely small contributions to rare decays, however in the presence 
of heavy KK neutrinos this is no longer the case. For example in the case of the original 
neutrino mass model of [2] the large splittings of the heavy neutrinos, together with their 
unsuppressed couplings to the SM fields yield a large loop- induced /x — > rate. In fact, 
Kitano [19] showed that a bound of ttikk > 25 TeV applies in this case. In our case however 
there is a generic softening of this bound, due to the fact that the SM fermions are now 
localized close to the UV brane, and therefore the charged current interactions with the KK 
neutrinos will be suppressed. This will happen generically in any model with bulk fermions. 
For the particular case of the A4 model the situation is even better: since the second Yukawa 
coupling of the charged leptons involving the right handed neutrinos is uniform due to the A 4 
representations, the interactions with all higgses (neutral or charged) will be diagonalized 
in the same basis where the masses are diagonalized. So one only needs to consider the 
exchange of charged W's and their KK towers, together with KK neutrinos. In fact, just as 
in the SM the loop induced contribution here will be finite. The reason for that is that the 
allowed additional brane localized counter term is of the form LHtfi' u^crF^, and this will 
be diagonalized once the charged lepton mass matrix is diagonal. This is again a specific 
property of the structure of the A 4 invariants. 

Next we will give a rough estimate of the KK mass bound from these processes. We will 
focus on loop induced /1 — > decays via to exchange of a W-bosons and KK neutrinos. The 
branching fraction from the exchange of a heavy neutrino and a W was calculated by Cheng 
and Li [26] and is given by (assuming the coupling to the W given by the usual SM gauge 
coupling g): 

3a 



B(fi -> ej) 



57T 



2 

m i 

Ml 



(5.1) 



where the sum over i indicates the sum over a generation of KK fermions, U is the PMNS 
mixing matrix between the SM charged leptons and a generation of KK neutrinos, and the 
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function F is given by 



F(z) = — 1 , , (10- 43* + 78* 2 - 49z 3 + 18z 3 log z + 4z 4 ). 
6(1 — zp 



(5.2) 



First we specify to the case of the exchange of a SM W and a KK mode neutrino. For 
z 3> 1 the approximate expression is F(z) ~ | + 3^p. Also, the coupling between the zero 
mode SM fermions, a KK neutrino and the zero mode W is suppressed at least by one factor 
of J'l, so there is an overall /| appearing in the rate. This is the main difference compared to 
the analysis of [19]: there all SM fermions were localized on the TeV brane, so the interaction 
with a KK neutrino was unsuppressed. The leading term drops out due to the unitarity of 
U, and so we are left with the approximation 

Of \ / M 54 " m W $ m KK i 2 m KK (t . Q x 

B(fi -> e 7 < fl f- — log , 5.3 

7T m\ K m z KK m w 

where Sttikk is the characteristic splitting among the KK modes in a given generation, given 
by « ^a-. We find, that even for Yukawa coupling close to the perturbative limit 

y ~ 3, and cl close to the composite case cl = 0.5 the branching ratio is two orders of 
magnitude below the experimental bound of 10~ n for a KK mass scale of 3 TeV. 

A slightly bigger contribution is obtained using the diagram where in addition to the 
KK neutrino one exchanges a KK W. In this case the gauge coupling could be as large 
as gfi^Jlog R' /R. One needs to take the function F(z) at values z ~ 1 for which it is 
approximately given by F(z ~ 1) ~ j| + J^(l — z). The universal piece drops out again due 
to the unitarity of U and we get an upper bound to the resulting branching fraction of order 

*(^7)</^^^Flog 2 f, (5.4) 
8007rm^ K m z KK R 

which again numerically is smaller than 10~ 13 for a 3 TeV KK mass. Clearly since the 
numerical contributions turn out to be not that far from experimentally interesting region 
it would be very interesting to perform a more detailed calculation of the // — > branching 
fraction in this model, including complete sums over KK towers (and also in general RS 
models with bulk fermions and Majorana neutrinos). 



6 Conclusions 

Warped extra dimensions provide a successful framework for flavor models: hierarchies in 
the masses and the mixing angles are naturally generated. Since neutrinos do not show 
hierarchies in the mixing angles, and only a mild hierarchy in the mass spectrum, one should 
introduce additional symmetries in the lepton sector. In this paper we have augmented the 
lepton sector of the RS model with the most successful and economical global symmetry used 
for neutrino mass models, the discrete non-abelian group A 4 . With appropriate assignments 
of the A 4 representations we can naturally achieve a successful lepton mixing pattern, while 
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the charged lepton mass hierarchy is implemented as usual in RS models via wave function 
overlap. The A4 symmetry also eliminates all tree- level sources of LFV in the neutral current 
sector, and so significantly lowers the bound on the KK mass scale. LFV appears only 
through charged current loops, and we estimated that the rate of /x — > is safely below the 
current experimental bound. So the most significant bounds on this model come from the 
EWP measurements, and as usual with KK mass scales of order 3 TeV these corrections will 
be also safely within the experimental bounds. 
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A Summary of A4 group theory 

A 4 is an 5*0(3) subgroup which leaves the tetrahedron invariant. It has 12 elements, two 
generators (S, T) connecting all of them, and four irreducible representations: one three- 
dimensional (3) and three one-dimensional (1, 1' and 1", with (1')* = 1"). Their products 
decompose as: 

V x 1" = 1, 1x3 = 3 

l'xl' = l", l'x3 = 3 ( , u 
1" x 1" = 1', 1" x 3 = 3 { ' 

3 X x 3 y = 3i + 3 2 + 1 + 1' + 1" 

where for the last line, with 3 X = (xi,x 2 ,x 3 ), 3 y = (2/1 ,2/2,2/3) and working in a basis where 
the three-dimensional representation of S is diagonal: 





AO \ 




(0 


1 




s = 


0-1 , 


T = 












\0 -lj 




V 








one has: 



1 = xm + x 2 y 2 + x 3 y 3 (A. 3) 

1' = Ziyi + uj 2 x 2 i/2 + ujx 3 y 3 (A.4) 

1" = xiyi + ujx 2 V2 + u 2 x 3 y 3 (A. 5) 

31 = (3:22/3,3732/1,^13^2) (A. 6) 

3 2 = (x 3 y2,x 1 y 3 ,x 2 yi) (A. 7) 
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with uj = e 27 ™/ 3 the cubic root of unity, satisfying: 

1 + uj + uj 2 = 0, uj* = uj 2 . (A.8) 

Note also that one has: 3 x 1' = 3 ~ u(xi,ujx 2 ,uj 2 x 3 ), where u ~ 1'. The same holds for 
3 x 1" with uj — > uj 2 . 

B 6*i3 = and 623 = 7r/4 at any order on IR brane 

We explicitly show in this section that the higher dimensional operators on the IR brane 
only affect 6w We recall that the Dirac mass matrix on the IR brane in presence of higher 
orders is of the form: 

(a P l\ 
7 a P] (B.9) 
P 7 ol) 

while at lowest order the Majorana mass matrix on the UV is: 

(a 0\ 

M U M = a d . (B.10) 
\0 d a) 

After the seesaw the Majorana mass matrix for the left-handed neutrinos (in the basis of 
diagonal charged leptons) shows the following pattern: 



Mu= c g f . (B.ll) 




Note that even for real input parameters this matrix has complex entries due to uj — e 2m / 3 
factors that do not cancel out when the left-handed doublet is rotated with V. If c and g 
were real, the Majorana mass matrix would be diagonalized with #13 = and #23 = 7r/4. 
In general it is a 3 x 3 complex symmetric matrix which thus contains 12 independent 
parameters. They are the 3 real eigenvalues, the 3 mixing angles and 6 phases. Moreover 
one can always redefine the neutrino fields to absorb 3 of them, leaving only 2 Majorana 
phases and a CKM one. Thus with the redefinition z/j — > e^z/j the mass matrix becomes: 

(^g2«<j!>l | c |g*(</ , c+<i!>l+02) | c |g*(</>l+02-0c)\ 

| c | e *Oc+</>i+<fe) |g| e i(2</> 2 +<j!> 9 ) Jgi(<fe+<fe) I _ (B.12) 

| c | e *0l+<?>2-</> c ) J e «(</>2+03) \g\ e i(2<t>3-(t>g) J 

It is now not difficult to see that this matrix can be made real with 0i = 0, 02 = —<p3 and 
03 = 0c provided that the relation 20 c = 9 holds. Although the expressions of these phases 
in terms of the input parameters are quite cumbersome, we checked that the latter relation 
is satisfied in our model. Therefore we conclude that the most general higher dimensional 
corrections on the IR brane can only modify 9 12 from its HPS value. 



18 



C Review of Gauge Breaking via BCs 



We shortly review here how the bulk SU(2) l X SU(2) R x U(l) b-l gauge symmetry is broken 
in this setup. The main goal is to define how to get the gauge boson profiles as we will need 
them to compute the W and Z couplings to the standard model fermions. We note A L , A R 
and A x as well as gm, g§R and g$ the gauge fields and coupling constants associated with 
SU(2)l, SU(2) r and U{1)b-l respectively. On the UV brane, the gauge symmetry breaks 
down to the SM gauge group SU(2) L x U(1)Y from Dirichlet BCQfor the SU(2) R fields: 

r A** = d M At* = 

[d z (g 5R A x + frA**) = 

where A ± = (A 1 TiA 2 )/V2. The chiral Slip's are broken to the vectorial subgroup on 
the IR brane by the finite VEV of a Higgs bidoublet (h) = diag(t>#, v H )/2, and the resulting 
(mixed) BCs are: 



r d z (g bL A L ^ - g bR Af a ) + V(g 5L A L ^ - g 5R A R ' a ) = 
z = R':{ d z (g 5R A L ^ + g 5L A^ a ) = (C.14) 

I d*A* = 

where V = (R' / R)(g§ L + g1 R )v 2 H /'i. The fifth components of the gauge fields have opposite 
BC. Obviously, these BCs allow for a massless (flat) mode which is nothing else but the 
photon field associated with the unbroken £7(1) of the (compactified) effective theory, and 
the KK decomposition is of the form: 

A^{x,z) = a L > R (z,m w )W±(x) + --- (C.15) 
A L ^ 3 (x,z) = J^a ^(x) + a L > R ' 3 (z,m z )Z fM (x) + --- (C.16) 

95L,R 

A*(x,z) = a 7(x) + a x (z,m z )Z ll (x) H (C.17) 

where the • • • stand for heavier KK resonances, and the wave functions are given by a(z, m) = 
z(AJx(mz) + BYiimz)). 

Yet remains the definition of the 5D gauge couplings. For this we have to match the latter on 
the 4D SM couplings. The fact that there is no SU(2) l x U(1)y symmetry in the effective 4D 
action makes the definition of the SM couplings somehow arbitrary. As matching conditions, 
we choose to recover the measured values of mw, mz and the electric charge e. After having 
fixed the values of R, R' and the ratio of the left/right gauge couplings, r = g^L/g^R, three 
parameters remain to be determined by the matching procedure, namely g$L, g$ and Vh- To 
fit them we proceed as follows. First, we relate g§ and g$i by imposing that the (canonically 



1 This can be thought as being the result of a Higgs mechanism in the limit where the localized scalar is 
decoupled. 
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normalized) photon couples to the electric charge Q = T 3 L + T 3R — Q B _ L /2. Given the KK 
decomposition above we get: 

1 (\ + r 2 1 \ 

+ ^)R\og(R'/R). (C.18) 



- 2 V all <%/ 

Then, the charged boson quantization equation fixes the product q^lVh as a function of 
my/- And plugging back these two relations into the neutral boson quantization equation 
allows to find the Higgs VEV as a function of mz, fnw an d e. Finally to fully determine the 
wave-functions we need to make the W and the Z are canonical fields by imposing: 

J' R dz(^j[a L ^(z,m w ) 2 + a R ^(z,m w ) 2 ] = 1 (C.19) 
j' R dz(^j [a L > 3 (z,m z ) 2 + a R > 3 {z,m z ) 2 + a x (z,m z ) 2 ] = 1. (C.20) 
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